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ABSTRACT 

We present a simple and efficient anisotropic generalization of the semi-isotropic (two- 
integral) axisymmetric Jeans formalism which is used to model the stellar kinematics of 
galaxies. The following is assumed: (i) a constant mass-to-light ratio M/L and (ii) a velocity 
ellipsoid that is aligned with cylindrical coordinates (R, z) and characterized by the classic 

anisotropy parameter Z = 1 — v^/v^. Our simple models are fit to SAURON integral-field 
observations of the stellar kinematics for a set of fast-rotator early-type galaxies. With only 
two free parameters (/3 2 and the inclination) the models generally provide remarkably good 
descriptions of the shape of the first (V) and second (V Ims = \/V 2 + a 2 ) velocity moments, 
once a detailed description of the surface brightness is given. This is consistent with previ- 
ous findings on the dynamical structure of these objects. With the observationally-motivated 
assumption that (3 Z > 0, the method is able to recover the inclination. The technique can be 
used to determine the dynamical mass-to-light ratios and angular momenta of early-type fast- 
rotators and spiral galaxies, especially when the quality of the data does not justify more 
sophisticated modeling approaches. This formalism allows for the inclusion of dark matter, 
supermassive black holes, spatially varying anisotropy, and multiple kinematic components. 
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- galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

According to the theory that best reproduces the observations, 
galaxy formation proceeds in a hierarchical fashion, driven by grav- 
ity, i n a Universe domin ated by dark matter of unknown nature 
(e.g. lSpringelit~aI].l2005h . The hierarchy of merging is accompa- 
nied by changes in galaxy structure and morphology. In particular 
early-type galaxies (Es and SOs) are thought to form by the gas-rich 
merging of spiral galaxies or by gas star vation of spirals, f ollowed 
by subsequent collisionless mergers (e.g. lFaber et all2007t) . 

Three key global parameters can be used to characterize galax- 
ies structure while studying this sequence of merging of galaxies 
and dark matter: (i) The angular momentum, which varies dur- 
ing mergers and increases with the amount of gas dissipation, (ii) 
the stellar population, which records the history of star formation 
events during the gas-rich mergers, and (iii) the mass-to-light ra- 
tio, which is affected by both the population and by the dark-matter 
fraction. 

The large majority of the galaxies in the Universe are to 
first order axisymmetric (except for bars) an d posses disks. 
This includes fast-rotator early-type galaxies ( Emsell em et alj 
2007; Cappella ri~et~all 120071) and spiral galaxies. Both the fast- 
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rotator early-typ es (Gerhard et al j|2001: iRu sin & KochaneM 2005 1 ; 
Cappellari et al. 2006; Koopmans et al. 2006; Thomas et al. 2007) 
and th e spiral galaxi e s |Persic et alj Il99rj; IPalunas & Wi lliams 
l200d : iBell & de Jon3 l200ll ; iKassin et alj I2OO6I) appear dom- 
inated by the stellar matter inside one half-light radius. 
Observations suggest that they have a dynamical structure 
characterized by a flattening of the velocity ellip soid in the z 
direct io n para ll el to the gala x y sym m etry axis (c ierssen et al. 



19971. l2000l: IShapiro et al.l 120031: ICappellari et al.1 
Noordermeer. Merrifield. & Arag6n-Salamancall2008h . 
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The goal of this paper is to include the knowledge of the struc- 
ture of the fast-rotator and spiral galaxies, into a realistic but sim- 
ple dynamical modeling method, which can be applied to the mea- 
surement of both the mass-to-light ratio and the amount of rotation 
(for which the inclination is needed) in the central regions of these 
galaxies, while also allowing for the inclusion of dark matter and 
the study of multiple kinematical components. The success of the 
adopted model's assumptions in describing the kinematics of real 
galaxies is verified against integral- field observations of the stellar 
kinematics obtained with SAURON dBacon et aill200lh . 

In Section[2]we briefly review the theory and past applications 
of the Jeans equations and of the shape of the velocity ellipsoid in 
galaxies. In Section [3] we describe our new anisotropic Jeans for- 
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malism, which we apply and test in Section [4] Finally our results 
are summarized in Section[5] 



2 SOLVING THE AXISYMMETRIC JEANS EQUATIONS 

2.1 The collisionless Boltzmann equation 

The positions x and velocities v of a large system of stars can 
be described by the distribution function (DF) /(x, v). When the 
system is in a steady state under the gravitational influence of a 
smooth potential $, the DF must satisfy the fundamental equation 
of st ellar dynamics, the st eady-state collisionless Boltzmann equa- 
tion jBinnev & Tremainell 19871 . hereafter BT; equation [4-13b]) 



I df _ 9£ df_\ = Q 
^ » ' dxi dxi dvi J 



(1) 



Given that / is a function of six variables, an infinite family of so- 
lutions satisfies equation (0. Additional assumptions and simplifi- 
cations are required for a practical application of the equation. One 
classic way of constraining the problem consists of drastically re- 
ducing it, from that of recovering the DF, to that of studying only 
the velocity moments of the DF. This approach leads to the Jeans 
equations, which are discussed in the next section. 



2.2 The Jeans equations 

2. 2. 1 Summary of derivation 

If we rewrite equation (T) in standard cylindrical coordinates 
(R, z, <fj) and we make the important assumption of axial symmetry 
(d§/d<j) = df/dct> = 0) we obtain (cf. BT equation [4-17]) 

df, df fvl 3*\ 3/ 9$ df vrv* df 
vrt^+v z - — h | — - | — — — 7^— = 0.(2) 



dR dz \ R dR I dvR dz dv z R dv$ 

Multiplying equation l|2} respectively by Vr and by v z , and int egrat- 
ing o ver all velocities, we obtain the twrj;] Jeans equations ( Jeans 
1922; BT equation [4-29a,c]) 
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where we use the notation 



uv k Vj = / v k Vjf d v. 



(3) 
(4) 

(5) 



These equations are still quite general, as they derive from the 
steady-state Boltzmann equation (TJ with the only assumption of 
axisymmetry. They do not require self consistency (a potential $ 
generated by the luminosity density and they make no assump- 
tion on the DF. However, even if one assumes $ to be known (it 
may be derived from the observed v via the Poisson equation), the 
two equations ([3j and © are still a function of the four unknown 
Vr, v z , v\ and vrv z and do not uniquely specify a solution. 



2.2.2 Closing the axisymmetric Jeans equations 

Given that the axisymmetric Jeans equations relate four functions 
of the meridional plane (R, z) coordinates, one needs to specify at 
least two functions of (R, z) for a unique solution. A natural way to 
provide a unique solution for the Jeans equations consists of spec- 
ifying the shape and orientation of the intersection of the velocity 
ellipsoid everywhere in the meridional plane. In fact the three com- 
ponents v R , v% and vrv z , of the velocity dispersion tensor, uniquely 
describe the equation of the velocity dispersion ellipse, whose shape 
can be derived by diagonalizing the tensor. The velocity dispersion 
ellipse will have the major axis at an angle 6 with respect to the 
equatorial plane 

2 v R v z 



tan 26» 



and an axial ratio < q < 1 given by 



Or+v* ~ V ( v r ~ v z) 2 +^v R v z 



(6) 



(7) 



+ v 2 z + J {v R - vl) 2 + 4 v R v z 



The specification of the ellipse geometry (8 and q) is then sufficient 
to write two of the three variables as a function of the remaining 
one. 

The most common choice which is made consists of assum- 
ing a circular velocity-ellipsoid in the meridional plane (some- 
times called semi-isotropy condition). This assumption implies 
v R — v z and vrv z = 0, and is sufficient to 'close' the set of 
equations to provide a unique solution for the remaining variables 
v z and v 2 . Other common options are discussed in Section l2l4l 

A well studied special case of semi-isotropic system is one in 
which the 3-integral DF which generally characterizes a stationary 
system, is assumed to depend only on the two classical integrals of 
motion / = f(E, L z ), where E is the potential energy, and L z 
is the angular momentum with respect to the symmetry z-axis. For 
this reason the semi-isotropic Jeans models are often called two- 
integral Jeans models. 

In this paper we are investigating whether it is possible to make 
an alternative assumption on the shape of the velocity ellipsoid, 
which retains the simplicity of the widely used semi-isotropic as- 
sumption, while also providing a better description of real galaxies 
and still leading to an efficient solution of the Jeans equations. 

2.2.3 Past applications of the semi-isotropic Jeans equations 

Due to its simplicity, the semi-isotropy assumption has proven 
remarkably useful in a large variety of applications of the 
Jeans equati ons. It was used to quant i fy the amount of rotation 
in ga l axies dNagai & Miyamoto] 1 19761 : ISatohl Il98d ; iBinnev et al] 
1990; van der Marel et al.l Il990h , to s e arch for hidden disks in 
ellipt ic al galaxie s dRix & W hite 1992; Cinzan o~& van der Marell 
Il994l; lRixetal.1 Il999l). t o measure t heir mass-to-light ra- 
tio ( Ivan der Marell Il99ll ; Istafler et al.l 1 19991 : ICappellari etal] 



2006; Cortes, Kenney, & Hardy 2008) and dark matter profiles 
i Carollo et alj 19951) , to study the connection b etween line-strength 
in galaxies and the local escape velocities jDavies et al.l 1 19931 ; 
Carollo & Danziger 1994), to s t udy t h e scatter in the Funda- 
mental Plane (van Albada et al.l 1 19951 : lLanzoni & Ciottil 120031 : 



1 All terms in the third Jeans equation, involving a multiplication by v^, 
are zero in axisymmetry. 



2 In the original paper by Jeans] j 19221) this was not an assumption. Stellar 
orbits were thought to conserve only two isolating integrals of motion, in 
which case the DF naturally possess that special semi-isotropic form. 
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IRiciputi et alj|2005l) . to measur e the deviations of the gas kine- 
matic s from circular velocit ies l Bertola et alj 19951 : ICorsini et alj 
1 19991 : ICinzano et all 1 19991 : lYoung et alj l2008h. and to esti- 
mate the masses of supermassive black holes (Magorrian et aL 
19981: Ivan der Marel et all Il998l : ICretton & van den Boschl 1 19991 : 



Josep h et alj2001 ). Although th ese models have been largely super- 



seded by [SchwarzschikJ d 19791) orbit-superposition method, when 
good kinematic data are available and the maximum general- 
ity is required, they are becoming useful to study the mass-to- 
light ratios and rotation of galaxies at high-redshift, where the 
data quality still does not justify mo re sophisticated approaches 
dvan der Marel & van Dokkumll2007allbT: ?). 



2.3 Shape and orientation of the velocity ellipsoid in galaxies 

To explore the possibility of making a simple but sufficiently realis- 
tic assumption on the shape and orientation of the velocity ellipsoid 
in the meridional plane of axisymmetric galaxies, we need to un- 
derstand what that shape is expected to be. 




0.2 0.4 0.6 0.! 

R 



2.3.1 Theory 

A qualitative insight into the orientation of the velocity ellipse in 
real galaxies can start from the an alysis of the special case of sep- 
arable potentials dde Zeeuwl[l985l) . In an axisymmetric oblate sep- 
arable potential the equations of motion for the stellar obits can be 
separated in a prolate spheroidal coordinates system (A, <j), v) which 
also defines the boundaries of the orbits. In other words the orbital 
motions can be written as the linear combination of two indepen- 
dent oscillations in A and v, plus a non-uniform rotation around the 
symmetry axis. For this reason the velocity ellipse is a ligned with 
the coordinate system and the cro ss term V\v v vanishes dEddingtorj 
Il915l : Ide Zeeuw & Hunter! 1 199Ch . The prolate spheroidal coordi- 
nates are characterized by the fact that they tend to be aligned with 
the cylindrical coordinates (R, z, (f>) at small radii and with the po- 
lar ones (r, 0, (f>) at large radii. Separable potentials are character- 
ized by a central constant-density region. The cylindrical alignment 
happens in that region, as the stars there tend to move like an har- 
monic oscillator. For these reasons one expects the velocity disper- 
sion ellipse to be cylindrically aligned at small radii and radially 
aligned at large radii. 

The gravitational potential of real galaxies is not of the separa- 
ble form. Separable potentials are in fact necessarily characterized 
by smooth analytic centers, while real galaxi es possess central sin- 
gularities due to c usped density profiles (e.g. lFerr arese et a l .11 19941: 
Lauer et al.lll995l) and supe r massive black holes dMagorrian et alj 
19981 : iGebhardt et alj I2OO0I : iFerrarese & Merrittl 120001) . Still, nu- 
merical integrations of orbits in non-separable axisymmetric poten- 
tials show that most of them are still bounded by cur ves, which 
qualitatively resemble the spheroidal coordinates ( e.g. iRichstonel 
1 1 9821 ; iDehnen & Gerhard] 1 19931 : ICretton et alj[l999l) . For this rea- 
son one can expect the velocity dispersion ellipse to behave in real 
galaxies in a way that is qualitatively similar to the separable case. 
Numerical calculations confirm this fact in non-separable triax- 
ial an d axisymmetric potentials (lMerrittlll980l : IDehnen & Gerhard! 
1993). In the limit of a point mass, the velocity ellipsoid has to be 
spherically aligned for symmetry. This suggests that at small radii, 
where the supermassive black hole or the stellar cusp dominates, 
the velocity ellipsoid in real galaxies should be more spherically 
oriented than in the separable case. 



Figure 1. Qualitative description of the shape and orientation of the velocity 
ellipsoid in the meridio nal plane (R, z) of the fast-rotator early-type galax- 
ies (derived from Cappellari et al. 2007). The ellipses show the intersection 
of the velocity ellipsoid with the (vr , v z ) plane. The solid lines show a 
representative prolate spheroidal system of coordinates. The dashed line is 
a representative iso-density contour. The ellipsoids are aligned in spheroidal 
coordinates, but the axial ratio of the ellipses varies with the polar angle in 
such a way that the shape and orientation of the ellipsoids are similar on 
both the equatorial plane and the symmetry axis. 

2.3.2 Observations 

The DF of a stationar y system i s a function of the three separable 
integrals of motion djeandfl9~15l) . In general it cannot be recovered 
from real galaxies without at least another three-dimensional ob- 
servable quantity. This quantity is now being provi ded by integral- 
field observations of the stellar kinematics (e.g. lEmsellem et alj 
2004), which allow the stellar line-of-sight velocity-distribution to 
be measured at every position of the galaxy image on the sky. 
We used these observations, in combination with orbit-based three- 
integral axisymmetric models, to measure the shape and orienta- 
tion of the velocity ellipsoid at every position within the meridional 
plane (within one half-light radius R a ), for a sample of 25 early- 
type galaxies (Cappe ilari et al.ll2007l) . Fig. Q] qualitatively summa- 
rizes the main findings of that paper regarding the shape of the ve- 
locity dispersion ellipsoid for the fast-rotator galaxies: (i) The ellip- 
soids qualitatively behave like in separable potentials and already 
near IR C they are essentially spherically oriented; (ii) the axial ra- 
tio of the ellipsoids varies gradually as a function of the polar angle, 
in such a way that the ellipsoid has nearly the same shape on both 
the equatorial plane (where the density is the highest) and the sym- 
metry axis. The net effect is that to first order the global anisotropy 
of the galaxies can be described as a simple flattening of the veloc- 
ity ellipsoid in the z-direction (w| < Vjj). 



2.4 Choice of the coordinate system for the Jeans solution 

One can think of three physically motivated choices of the coordi- 
nates system in which to align the velocity ellipsoid, for a unique 
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solution of the axisymmetric Jeans equations. In this section we ex- 
amine advantages and problems of each one in turn. 



2.4.1 Prolate spheroidal coordinates 

The observed behavior of the ellipsoid could be approximated 
by solving the Jeans equations in p rolate spheroidal coordinates 



in a generic axisymmetric potent ial (Dejonghe & de Zeeu wll 19881 : 
lEvans & Lvnden Bell|[l99"ll : lArnoldll995l : Ivan de Ven et alj|2003h . 



To qualitatively reproduce the behavior of Fig. Q] the simplest real- 
istic model would need at least four free parameters. Two param- 
eters are required to define the coordinate system. A third param- 
eter could define the shape of the velocity ellipse on the equato- 
rial plane and another parameter would describe the angular varia- 
tion. A problem of this approach is that the general solution, to de- 
rive the observables for realistic galaxy potentials, requires at least 
a computationally-expensive triple0 numerical quadrature. The de- 
rived equations are also rather cumbersome and not justified by the 
fact that these models can at best provide a qualitative description of 
real galaxies. For these reasons those solutions have currently only 
been applied to a handful of analytic potentials, and no application 
to real galaxies exist. 

2.4.2 Spherical coordinates 

Given that the velocity ellipsoids in Fig. [T] are nearly spheri- 
cally aligned, a much simpler choice would consist of solving the 
anisotropic Jeans equation i n spherical coordinates. A solution for 
this case w as presented b y [Bacon et aU dl983|) and applied to real 
galaxies bv lBacor] Jl985h and lFillmord il986l) . However even un- 
der this spherical assumption the solution still has the same form as 
in the spheroidal coordinates. The computation still involves a triple 
numerical quadrature. This is in fact a generic feature of the axisym- 
metric Jeans solution: even in the simple case in which the density is 
assumed to be stratified on similar oblate spheroids, one quadrature 
is needed to obtain the potential, a second quadrature provides the 
intrinsic velocity's second moments and a third quadrature finally 
gives the projected observables. In the spherically-oriented case the 
coordinates system has no free parameters and the simplest realis- 
tic model to reproduce Fig.Q]could have two parameters (anisotropy 
and its angular variation). 



2.4.3 Cylindrical coordinates 

A more radical alternative is to assume the velocity ellipse is 
oriented with the cylindrical c oordinate system. This option was 
tested by e.g. Fillmore] Jl986h . however until now, for the rea- 
sons described in Section 12.4.11 it was not considered a sensi- 
ble choice to describe the shape of the velocity ellipsoid in real 
galaxies. O ur new results, based o n the SAURON integral-field ob- 
servations ( Cap pellari et al. I l2007t) , make this choice worth explor- 
ing again. This option cannot provide a formally accurate rep- 
resentation of Fig. Q] however it is accurate near the equatorial 
plane, where the density is at its maximum, and near the minor 
axis. Models with cylindrically-oriented velocity ellipsoid provide 
a good qualitative description of the empirical observation that 
the global anisotropy in fast-rotator galaxies is best characterized 
as a flattening of the velocity ellipsoid in the vertical ^-direction 

3 For typical accuracies, every additional numerical quadrature generally 
increases the computation time by about two orders of magnitude. 



( Cappellari et al. I l200l their Fig. 2). The near-cylindrical orienta- 
tion of the velocity ellipsoid may be due to the presence of disks, 
where this orientation appears natural. For this reason the cylin- 
drical orientation i s certainly appropria t e to describe the dynamics 
of spiral galaxies jGerssen et al.lll997l. |2000|: IShapiro et alJliool : 



iNoordermeer. Merrifield. & Arag6n-Salamancal f2008). 

In the next section we show that the assumption of 
a cylindrically-aligned velocity-dispersion ellipsoid, combined 
with the powerful Mu lti-Gaussian Expansion (MGE) method of 
lEmsellem et alj Jl994) . can generate simple solutions that well re- 
produce the integral-field kinematics of real galaxies, and also allow 
for variable mass-to-light ratios (e.g. dark matter), spatially varying 
anisotropy and multiple kinematical components. All this while still 
requiring a single numerical quadrature to predict the observables 
on the sky plane (or any other projection). 



3 ANISOTROPIC JEANS SOLUTIONS 
3.1 Axisymmetric case 

3.1.1 General solution 

We start from the general axisymmetric Jeans equations {3]> and © 
and we make the following two assumptions: (i) the velocity ellip- 
soid is aligned with the cylindrical coordinate system (R, z, (f>) and 
(ii) the anisotropy is constant and quantified by v\ = b w|. In this 
case the Jeans equations reduce to 



bvvl 



R 



+ 



d(b vvj) 
dR 

dz 



'dR 
' dz' 



(8) 



(9) 



which corresponds to the semi-isotropic case (two-integral) when 
6 = 1. With the boundary condition vv% = as z — > 00 the 
solution reads 



uv%(R,z) = 



dz 

R ^R~ +VV * 



p d ® 
d~R' 



(10) 



(11) 



A general caveat regarding the Jeans equations is that the existence 
of a solution does not guarantee the existence of a corresponding 
physical positive DF. This can only be verified using different tech- 
niques. 

3.1.2 Summary of MGE formalism 

To derive solutions for the Jeans equations we make an explicit 
choice for the parametrization of the stellar density and the total 
density (which can include dark matter and a cen tral black hole). 
We ad opt for both the MGE parametrization of lEmsellem et al] 
Jl994l) due to its flexibility in accurately reproducing the surface- 
brightness of real galaxies and the availabil ity of robust routine jf| to 
fit the galaxy photometry dCappellaril2002l) . If the x-axis is aligned 
with the photometric major axis, the surface brightness E at the lo- 
cation (a/, y') on the plane of the sky can be written as 



T,{x ,y' 



exp 



1°l 



12 

12 . y 
x + — 



(12) 



4 Available from http://www-astro.physics.ox.ac.uk/~mxc/idl/ 
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where N is the number of the adopted Gaussian components, hav- 
ing total luminosity Lk, observed axial ratio < q' k < 1 and dis- 
persion o-fc along the major axis. 

The deprojection of the surface brightness to obtain the intrin- 
sic luminosity density is not u nique unless the axisymmetric galaxy 
is see n edge-on (t = 90°) dRvbickil [l987l ; iKochanek & Rvbickil 
1996), and the degeneracy becomes e specially severe when the 
galaxy is seen at low inclinations dGer hard & Binney 1 19961 : 



Romanowskv & Kochaneklll997l : Ivan den Boschlll997| ; Magorrian 
1999h . The M GE method provides a simple possible choice for the 



deprojection dMonnet et al.l 1 19921) . One of the advantages of the 
MGE method is that on e can easily enforce the roundness of the 
model l lCapDellarill2002n . thus producing realistic densities, which 
look like real galaxies when projected at any angle. The method 
cannot eliminate the intrinsic degeneracy of the deprojection and 
this has to be considered when interpreting results of galaxies that 
are close to face-on. The deprojected MGE oblate axisymmetric lu- 
minous density v can be written as 



(V27T Ofc) 3 <?fc 



exp 



2*2 



R 2 + Z - 

4 



(13) 



where the individual components have the same luminosity L k and 
dispersion a k as in the projected case J 1 2b . and the intrinsic axial 
ratio of each Gaussian becomes 



5fc = 



(14) 



where i is the galaxy inclination (i == 90° being edge-on). 

The total density p can be generally described by a different 
set of M Gaussian components 



p(R, 



E 



Mi 



27T a j) 3 qj 



exp 



1 

-2^ 



„9 Z 

R + - 



(15) 



In the self-consistent case the Gaussians are the same as in equa- 
tion U3t and one has M = N, Oj = <Tk, Qj — Qk and Mj = T k L k , 
where T k is the mass-to-light ratio, which can be different for dif- 
ferent components. In the non-self-consistent case the density can 
be described with the sum of two sets of Gaussians: the first de- 
rived by deprojecting the surface brightness with equation H3\ , and 
the second e.g. obtained by fitting a (one-dimensional) MGE model 
to some adopted an alytic parametrization for the dark matter (e.g. 
lNavarroetal]|l996l) . 

The gravitationa l potential generated by the density of equa- 
tion dl5t is given by jEmsellem et alJI 19941 ) 




*(«,*) = - 



where G is the gravitational constant and with 



(16) 



exp 



R 2 + 



1-(1 



q 2 )u 2 



(17) 



A supermassive black hole can be modeled by adding a Keplerian 
potential 

GM. 



(R, z) = 



(18) 



VR 2 + z 2 

to equ ation H6\ . However, as pointed out by lEmsellem et alj 
(1994), an even simpler approach consists of modeling the black 
hole as a small Gaussian in equation l !16t . having mass Mj = M,, 



qj — 1 and 3<Tj < r nl - ln , where r m j n is the smallest distance from 
the black hole that one needs to accurately model (e.g. one could 
choose r m i n w cr ps f). 



3.1.3 MGE Jeans solution 

Now we apply the MGE formalism to the solution of the axisym- 
metric anisotropic Jeans equations of Section 13.1.11 Our deriva- 
tion is an extension of what was done in the semi-isotropic self- 
consistent case (bk_= 1 and Mj = TL k ) i n section 3. 4 of 
lEmsellem etai] ll994h . As already pointed out bv |jeanslfl922T) . his 
equations can be used to model the kinematics of different dynam- 
ical tracers, as long as they all move in the same potential (e.g. 
iRix & White] 1 19921 : ICinzano & van der Marellll994l) . To maintain 
generality we will then write the solution for the individual N lu- 
minous Gaussian components, which can then be assumed to have 
different anisotropy. This fact can be used e.g. to model anisotropy 
gradients, or to study the anisotropy of kinematical subcomponents 
in galaxies (e.g. bulge and disk). Substituting equations d 1 3 1 > and 
l !16t into equations j 10b and i ll IK the integral in z can be performed 
analytically and we obtain 



» = 47rG / E T-Cu 2 



lM ~ 2 " 2 *-q jP vHj{u)u> 



da 



W v l\k = b k [vv 2 ] k 



1 v-~* VkqjPojHj (m) u 2 2 



./n . -. ' 



3=1 



Cu 2 



VR'&u 



(19) 
(20) 

(21) 



1 -^VkqjPojHjiuju 2 , 2 2 2 



.In . . 1 



3=1 



Cu 2 



VR +b k a k q k )Au, 



where we defined v k — v k (R, z), poj = pj(0, 0) and 
C = 1 



2 

9i 



2 2 

°k ik 



v 



b k ql- [(l~b k )C + {l~q 2 )b k ] 



(22) 
(23) 



In all the equations of this paper the index k refers to the parameters, 
or the anisotropy, of the Gaussians describing the galaxy's luminos- 
ity density (equation 1131 ), while the index j refers to the parame- 
ters of the Gaussians describing the total mass (equation 1151 ). from 
which the potential is obtained. 

When b k is not the same for the individual luminous Gaus- 
sians, the total luminosity-weighted anisotropy at a certain spatial 
location (R, z) of an MG E model is given by the standard definition 
dBinnev & Mamonlll982l) . combined with equation H9i : 



/3 z (R,z) 



Y,k bk[w 2 ] k 



J2k bkVk 



(24) 



The last approximation comes from the fact that [v 2 ] k , being mostly 
a function of the total MGE potential, varies relatively little for the 
different Gaussians, while v k can be completely different and varies 
by many orders of magnitude for the various luminous MGE com- 
ponents. This allows the global anisotropy of an MGE model, at a 
certain spatial location in the meridional plane, to be approximately 
estimated from a simple luminosity- weighted sum of b k . A similar 
reasoning applies to the estimation of the total k parameter in in 
Section l3.1.5l and the total fi parameter in Section l3.2.2l 
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3.1.4 Line-of-sight integration of the velocity second moment 

The intrinsic quantities have to be integrated along the line-of-sight 
(LOS) to generate the observables which can be compared with the 
galaxy kinematics. For this we define a system of sky coordinates 
with the z axis along the LOS and the x' axis aligned with the 
galaxy projected major axis (see equation 1121 ), The galaxy coordi- 
nates [x, y, z) are related to the ones in the sky plane by 



(25) 
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sini 
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sin i 


cos i 
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where the z-axis coincides with the galaxy symmetry axis and the 
cylindrical radius is defied by R 2 
ond velocity moment along the line-of-sight v? B 
luminous Gaussian component, is then given byj 



x 2 + y 2 . The projected sec- 
for one 



los-l 



+ 



{ 



VVi Ife cos % 



yv 2 R ] k sin 2 <t> + [vv 2 ^ cos 2 <p \ sin 2 i\ 



dz , 



(26) 



where cos <j> = x/R, while the total observed second moment, for 
the whole MGE model, will be 



£«, 2 

los 



(27) 



After substitution of equations J 191 — d2 lb . the z integral can be writ- 
ten explicitly. Summing over all the N luminous Gaussian compo- 
nents we obtain the final expression 

r i N M 
4tt 3/2 G / Vak * P°i ^ 

fe = l j = l 

2 , 



o\q\ (cos 2 i + bk sin 2 i) + T)i 



(1 - Cv?) J (A + B cos 2 i) [1 
(A + B)y' 2 



1 



q 2 )u 2 



x exp 



-A 



where we defined vok 



A + B cos 2 i 
u k (0,0) and 



du, 



A = 



+ 



'it 



+ 



(1 



? 2 )« 4 



(1 



q 2 )u 2 ] 



(28) 



(29) 



(30) 



As expected e quatio n i2S\ reduce^] to equation (61) of 
lEmsellem et"ai] i 1994) when bk = 1 and Mj = TLk- As in 
the semi-isotropic case this formula is very quick to evaluate as it 
still requires a single numerical quadrature and involves no spe- 
cial functions. All this starting directly from a fit to the galaxy 
surface-brightness, without the need for numerical deprojection or 
PSF deconvolution, as required in other approaches. The formula 
has various possible applications, as it can be used to model realis- 
tic galaxies with variable anisotropy or multiple kinematic compo- 
nents, dark-matter, variable stellar M/L and supermassive black 



5 Completely analogue expressions can b e found for the proper m otions, 
using e.g. the formulas in appendix A of [Evans &deZeeuwH1994l) . We 
found that all three components of the projected proper motion dispersion 
tensor c an be written via single quadratures and using no special functions. 

6 As in Cappellari et al. (2006) we corrected the typo in the expression for 
B in equation (63) of Emsellem et al. ( 1994). 



holes (representing the point mass with a small gaussian as de- 
scribed in the last paragraph of Section [3.1.2t . 

Simple expressions, involving a single quadrature, can be de- 
rived for the second velocity moments also for the case of the Ke- 
plerian potential l l 1 8t of a black hole, as done in the semi-isotropic 
case in appendix A of Em sellem et alj d 1994b . The black hole mo- 
ments can then be quadratically co-added to the ones in equa- 
tion d28t to obtain the observed velocity moments for the galaxy. 
We found no advantage in speed or accuracy when performing a 
separate calculation for the black hole and the galaxy potential. For 
this reason we will then not give separate expressions for the Ke- 
plerian case. Using equation d28b with a black hole it is however 
important to use a quadrature routine which samples the integrand 
function at a sufficiently high number of initial points, to properly 
recognize the sharp peak in the integrand near u = 0. Here we us ed 
the vectorized adaptive quadrature algorithm o f lShampinel ( l2008l) . 

The second moments u 2 os provided by equation J28t are a good 
approximation for the observed quantity V 2 ms = V 2 + a 2 , where 

V is the stellar mean v elocity and a is the velocity dispersion. In 
Cappellari et we used realistic semianalytic dynamical 
models of galaxies and found that, to extract vf oB from the simu- 
lated data, one should use a single Gaussian LOSVD and adopt as 

V and a the mean velocity and dispersion of that Gaussian. Due to 
the sensitivity of the second moments to the uncertain wings of the 
LOSVD, this approach is preferable than trying to extract a more 
compl ex LOSVD, e.g. by fitting the Gauss-Herm ite parametriza- 
tion ivan der Marel & Franxll 19931 : lGerhardll993l) or a fully non- 
parametric LOSVD, and numerically integrate vf oB from that. 



3.1.5 Line-of-sight integration of the velocity first moment 
The projected first velocity moments v\ai = : *V are given by 



VV,f, cos 



isini dz' . 



(31) 



In this case, the two assumptions we made in Section [3.1.1l are not 
sufficient any more to provide a unique prediction and therefore 
additional assumptions are needed. The Jeans equations and ((9) 
in fact only give a prediction for and one has to decide how 
the second moments separate into the contribution of ordered and 
random motion, as defined by 



— 2 2 . 2 



(32) 



This need for extra assumptions on the tangential anisotropy is a 
fundamental limitation of the first-moments equations, and it is the 
reason why one should fit the more general equation J28 1 to V rms , 
and only subsequently fit any extra parameter of the first moment 
solution to V, instead of simultaneously fitting the Jeans solutions 
to both V and a. 

The first moment equations however are very useful to quan- 
tify the amount of rotation in galaxies and for this reason have been 
used in the past (Section [2.2.3t . One can think of two natural op- 
tions to parameterize the separation of random and ordered stream- 
ing rotation around the symmetry axis. The first option consists of 
assuming a constant anisotropy, for each Gaussian component, in 
the (vn,Vj,) coordinates, analogously to equation d!9t 



which would imply 



[IV, 



\l/2 

c k [v%]k 



(33) 



(34) 
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A practical limitation of this approach is that the mean velocity 
would depend in a non-linear way on the c k parameter. For this rea- 
son nearly all previous authors ado pted, in the se mi-isotropic case, 
a second approach, introduced by ISatohl {1980). This consists of 
defining a constant k which quantifies how much the model velocity 
field is scaled with respect to that of the isotropic rotator (|fc| = f). 
In this way the first moment can be computed only once, for the 

isotropic rotator, and is scaled linearly with k. 

W e adopt here the analogue, in our anisotropic case, of [Satoh 
(1980) approach. It has the added advantage of providing a direct 
measure of the amount of rotation, which is a more meaningful 
quantity than the tangential anisotropy. We define for each Gaus- 
sian component 



/ \l/2 

Mfc = K-k [\vl\k ~ [V%] h j 



(35) 



Here = when the fc-th Gaussian component is not rotating 
and = 1 when its velocity ellips oid is a circl e in the (vr, v<f,) 
plane. If bk = 1 then Kk reduces to ISatohl dl980h parameter and 
in this case \nk\ = 1 implies isotropy (the velocity ellipsoid is a 
sphere everywhere). An upper limit to \nk\ is set by the physical 
requirement that <r^ > (equation 1321 ). When summed over all 
luminous Gaussian components of the MGE modeQ 



3.2 Spherical case 

The cylindrically oriented assumption for the shape of the veloc- 
ity ellipsoid discussed in Section [37X1 is likely to be unrealistic for 
slow-rotator elliptical galaxies. As a class these objects are weakly 
triaxial, and i ndeed some of them show clear twists in their kine- 
matical axes ( Kormendy & Bender 1996; Emsell em et al. 1 120071 : 



ICappellari et alj2007l : lKrainovic et alj2008h . The most general and 
realisti c models for these objects are triaxial and orbit- or particle - 
based Ide Lorenzi et alj|2007l : lvan den Bosch. & van de Venl 2008). 
However also in the triaxial approximation a degeneracy in the 
recovery of the galaxy shape is still generally present so that no 
unique solution can be obtained. When one is only interested in 
global galaxy quantities, or to test the results of more general mod- 
els, it is still useful to construct simpler and approximate models. 
The isophotes of the slow rotators are in projection close to circular, 
especially in their central parts where the kinematics is generally 
obtained, and this implies they must be intrinsically not far from 
spherical. This suggests that one can use simple spherical models 
as a first order approximation to the dynamics of at least some of 
these objects. The spherical solution of the anisotropic Jeans equa- 
tions, in the MGE formalism, will be discussed in the following 
sections. 



W4, 2 = ^[i/i^ 2 ] fc , 
fe=i 

so equation ( 1351 ) impliefl 

N 

E 



fc = l 



1/2 



(36) 



(37) 



Substituting equation ([37) into equation <3U , and using equa- 
tions (|19M21b. we obtain the projected first velocity moment of the 
whole MGE model 



E 11103(2;', y) 



2V tiG x sini 

1 N M 



(38) 



EE 

fc = l j = l 



K k u k qjPQj'H.j(u) u V 
l-Cu 2 



du 



1/2 
dz'. 



When both bk = \kh\ = 1 and Mj — TL k , this equation re- 
duces to a self- c onsist ent isotropic rotator as in equation (59) of 
lEmsellem et alj dl994l) . A double quadrature seems unavoidable 
here, but when k is assumed to be constant for the whole MGE, 
this integral has to be evaluated only once with K k = 1, at the best 
fitting (i, /3 Z , T) parameters previously determined from a fit to the 
more general second moment equation d28t . and then vu^ can be 
linearly scaled by k to fit the data0 



7 As v^ 2 is denned as a quadratic difference of second moments it is not 
correct to write livZ = X^fcLi [ u7 "4>\k> as one would normally expect for a 
first moment. 

8 To allow for counter-rotating Gaussian components, namely to keep track 
of the sign of K k in equation and <38h one should multiply each 
term of the fc-summation by sgn(n k ) = Kk/\n^\ and then compute 
sgn(ui) l^l 1 / 2 , where w = [...] is the expression inside the big square 
brackets. 

9 When the alternative definition for the tangential anisotropy of equa- 
tion (33) is used instead of equation (35), the term T> in equation (38) be- 
comes [V+ (1 - cj.) b k qla 2 k /R 2 ]. 



3.2.1 General solution 

When the Boltzmann equation (JJ is written in spherical coordi- 
nates (r,9,<f>), by analogy with our derivation of the axisymmet- 
ric Jeans eq uations one can obtain t he Jeans equations in spherical 
symmetry (Bin nev & Mamonll 19821 ; equation [4-54] of BT) 



dr r ar 



(39) 



where v 2 = w? for symmetry and we defined j3 = 1 — v 2 /v 2 . The 
solution of this linear first-order differential equation with constant 
anisotropy /3 and the boundary co ndition vv 2 = as r — > 00 is 
given by (e.g. Ivan der Marelll 19941) 



vv 2 (r) 



-20 



[ 



Gr 



-10 



V[U) : U du 

du 

u{u)M{u) 



..2-20 



du, 



(40) 



considering that d$/dr = GM/r 2 . After projection along the line- 
of-sight z' , we obtain the observed second velocity moment (see 
equation [4-60] of BT) 



= 2G 



4 



vv 2 (r) 
\fr 2 



R 2 



: dr 



R 2 0) f°° u(u)M(u) 



v7 



R 2 



.2-20 



du 



(41) 



dr, 



where R is the projected radius, measured from the galaxy center. 
Integrating by parts the double integral can be reduced to a single 
quadrature, involving special functions 

3 T(r)is(r)M(r) _ 

R 



2G 



r 2-20 



-dr, 



(42) 



where 



Hr) 



R 



1-20 



/?(§ 



(3B W 1/3 

/3)r(/3-i) 



1 1 

2' 2 



B w f3 



1 1 

2' 2 



(43) 
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with w — (R/r) 2 , F is the Gamma func tion and B w is the in- 
comp lete Beta function (equation [6.6.1] of lAbramowitz & Stegunl 
1965). In the isotropic limit 



lim T(r) = Vr 2 ~R 2 , 



(44) 



and equation J42b reduces to the spherical isotropic form of equa- 
tion (29) of iTremaine et alj dl994f) . Some numerical implementa- 
tions provi de the incom plete Beta function only for positive argu- 
ments (e.g. IPress et al.M 1992r> . for negative values (/3 < 1 /2) one 
can use its analytic continua tion via Gauss's Hypergeome tric func- 
tion 2 Ji (equation [6.6.8] of lAbramowitz & Stegunll 19651) 



B w (a,b) = — 2-F1 (a, 1 — b; a + 1; w), 



(45) 



for which efficient routine^] exist dShaniie & Jianming 1996]). For 
(3 — ±1/2, the B w function is divergent, but for all practical pur- 
poses we found it is sufficient to perturb the {3 value by a negligible 
amount to avoid the singularity. 



3.2.2 MGE spherical Jeans solution 

We proceed as in Section |3.1.3| to derive an explicit solution for the 
spherical Jeans equation using the MGE parametrization for both 
the luminosity density and the total density. The expressions for the 
surface brightness, luminosity density and total density are as i n 
equation G2), GU and (B), with q 3 = q k = 1 dBendinelli|[l99lh . 
For a single Gaussian component they read 



E fc (ii) = 



2nol 



exp 



Li 



{V2^a k ) 3 

Mi 



exp 



2*2 



■ exp 



2a] 



(46) 



(47) 



(48) 



The mass of a Gaussian contained within the spherical radius r is 



M i {r) = M j 



erf 



s/2i 



■ exp 



2o 2 



(49) 



with erf(x) the error function (equation [7.1.1] of 
lAbramowitz & Stegunl 1 1965b . Using equation J42t and sum- 
ming over all Gaussian components with equation 127b . we obtain 
the projected second velocity moment for the whole MGE model 



2G 



E 



^k(r) v k (r) 

j,2 — 2/3 k 



dr, (50) 



where Vk(?) and Mj(r) are given by equation J47b and ( 1491 ) re- 
spectively, and J~k{r) is obtained by replacing the /3 parameter in 
equation d43 b with the anisotropy {3k of each luminous Gaussian 
component. We explicitly included the mass M, of a central super- 
massive black hole. As in the axisymmetric case, this formula in- 
volves a single numerical quadrature. It can be used to model nearly 
spherical objects with a variable anisotropy profile, variable stellar 
M/L, a supermassive black hole and dark matter. 

An approximate way to construct a model with a certain 
smooth radial profile of anisotropy /3(r), consists of defining [3k = 



f3(ok) for each MGE Gaussian with dispersion a k - This simple 
method works because in a spherical MGE each fc-th Gaussian pro- 
duces a sig nificant contributi on to the total luminosity density only 
near its Ok dCappellarill2002l ; see also equation 1241 ). 



3.3 Availability 

The Jeans Anisotropic MGE (JAM) package of IDLFl procedures, 
providing a reference implementation for the equations described 
in this section, together with other routines to evaluate auxiliary 
quantities, like the circular velocity, from MGE models, is available 
online from http://www-astro.physics.ox.ac.uk/~mxc/idl/ 

4 APPLICATIONS AND TESTS 



The formalism derived in Sections 13.1.41 and 13.1.51 generalizes to 
anisotropy the widely used semi-isotropic formalism for the so- 
lution of the Jeans equations while maintaining its simplicity. It 
was motivated by our model ing results on the orbi tal distribution 
of a sample of real galaxies ( Cappellari et al. 2007). However the 
method would be of limited usefulness if it did not describe the 
dynamics of fast-rotator early-type galaxies well. In this section 
we show that the simple assumptions we made provide a remark- 
able good description of the main features of the kinematic!* 2 ! 
dEmsellem et al.l2 004) of the fast rotator early-type gal axies, as ob- 
serve d with the SAURON integral-fi eld spectrograph ( Bacon et alj 
l200ll) as part of the SAURON survey dde Zeeuw et alj|2002h . 

4.1 Comparison with more general models 

4.1.1 Anisotropy comparison 

In this section we compare the anisotropy derived from the Jeans 
models, with the results for the glob al anisotropy obtained with the 
more general [Schwarzschil dl d 1979b method on th e same galaxies 
and from the same data in Cappellari et alj d2007I) . We use equa- 
tion j28t , with constant anisotropy 6 fc = 6 and constant mass-to- 
light ratio T for the whole MGE model, to predict the velocity 
second moments. Following standard practice we parameterize the 
anisotropy with the variable [3 Z = 1 — v\lv\ = 1 — 1/6. The 
problem becomes a function of two nonlinear variables (i,f3 z ) and 
of the scaling factor T. To further limit the number of free variables 
in the model, and to eliminate the degeneracy in the deprojection of 
the surface brightness, we select for the comparison the galaxies for 
which the photometry already constrains the inclination to be close 
to edge-on (i = 90°). For this we selec t from the MGE models 
in Table B1-B2 of Cappella ri et al.1 d2006i) . the fast-rotator galaxies 
for which the flattest Gaussian has q < 0.35, which forces a strict 
limit on the inclination i > 70°. We further exclude the galaxy 
NGC 4550, due to the presence of two counter-rotating stellar disks, 
which complicate the interpretation of the models, and NGC 4526, 
which has a strong dust disk affecting the observed stellar kinemat- 
ics. This selection leads to the four galaxies NGC 821, NGC 3377, 
NGC 4621 and NGC 5845. 

We evaluated the MGE model predictions given by equa- 
tion J28b . adopting i = 90 as in ICappellari et al.1 (2007), for a 
direct comparison with the anisotropy determinations. The model 
predictions were convolved with the PSF and integrated over the 



10 Available from http://jin.ece.uiuc.edu/routines/routines.html 



11 http://www.ittvis.com/idl/ 

12 Available from http://www.strw.leidenuniv.nl/sauron/ 
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NGC0B21 NGC3377 NGC4621 MGC5B45 
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Figure 2. Data-model comparison for the second velocity moments of edge-on galaxies. From left to right, the columns show the galaxies NGC 821, NGC 3377, 
NGC 4621 and NGC 5845. The top row visualizes the bi- symmetrized SAURON observations for V rms = VV 2 + a 2 . The subsequent rows show the model 
predictions for (" 2 os ) 1/ ' 2 as given by equation (28) for (3 Z = (semi-isotropic model; second row) and anisotropic models with f3 z = 0.1, 0.2 and 0.3 
respectively. The color levels are the same for the five panels in each column. NGC 4621 appears well reproduced by a weakly anisotropic model (3 Z 0.1. 
NGC 812 and NGC 5845 are well described by a model with /3 Z «i 0.2, while only the strongest anisotropy /3 Z 0.3 can qualitatively describe the observed 
shallow gradient of Vims along the projected minor axis of NGC 3377. These estimates are more accurately quantified in Fig. [5] 



pixels, before comparison with the observables, as described in Ap- 
pendix [A] For each galaxy we computed the model predictions for 
(^ios) 1 ^ 2 inside each Voronoi bin on the sky, at different anisotropy 
values p z . At every anisotropy the best fitting mass-to-light ratio T 
is obtained from the simple scaling relation T oc vf os as a linear 
least-squares fit. The best-fitting scaling factor for the model ve- 
locities is s — cos a |d|/|m|, where a is the cosine of the angle 
between the data d and model m vectors, which implies P^l 



Anisotropy NGC0821 



Anisotropy NGC3377 



T 



d ■ m 



m • m 



(51) 



where the vectors d and m have elements d n = [V r ms/AKma]n 
and m„ = [^L) 1/2 /AV lma ] n , with [V lns ]„ = y/V 2 + al 
the me asured values for the P Voronoi bins from Emsellem et al. 
(2004) and [AV rms ]„ the corresponding errors, while [(v 2 -,,;) 1 ^ 2 ]™ 
are the model predictions for T = 1. A qualitative comparison be- 
tween the observed velocity second moments and the model ones, 
for different values of the anisotropy parameter f3 z is shown in 
Fig.[2] 

To ease the visual comparison, given that the models are 
bi-symmetric by construction, the observations have been bi- 
sy mmetrized with respect to the kinematical major axis PA^m given 
in Cappella ri et alj d2007T) . The actual fits and \ 2 calculations are 
always performed using the original data and errors. It is clear that 
the isotropic models (J3 Z — 0) do not provide a good description 
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Figure 3. Best-fitting Jeans anisotropy. A\ 2 = X 2 ~ X m i n - describ- 
ing the agreement between the data and the models, is plotted against the 
anisotropy parameter fi z . The four panels show the galaxies NGC 821, 
NGC 3377, NGC 4621 and NGC 5845. The solid line with open circles rep- 
resent the Ax 2 of the Jeans models. The dashed horizontal line indicates 
the level A\ 2 = 9, which corresponds to the 3cr confidence level for one 
degree-of-freedom. The filled square with error bars shows the value of the 
global anisotropy measured using Schwarzschild models in Cappell ari et al l 



13 We corrected a typo in the corresponding equation (2) of Cappellari et al. 
j2006h . 
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for the observed second velocity moments of any of the four galax- 
ief^l the disagreement being strongest for NGC 3377 and weak- 
est for NGC 4621. Using the models presented in this paper we 
can now move away from the isotropic assumption and quantify 
that the model of NGC 4621 requires a small amount of anisotropy 
to appear qualitatively like the data: any /3 Z > 0.1 produces in fact 
two vertical 'lobes' along the minor axis, which are not observed. 
Both NGC 821 and NGC 5845 require a more significant anisotropy 
flz ~ 0.2 to reduce the amount of rotation (high [« 2 os ] 1//2 ) along 
the major axis to the level of the observations. A larger (3 Z ~ 0.3 is 
strongly excluded by the data, as the model immediately produces a 
significant elongation along the minor axis, which is not observed. 
Finally, in the case of NGC 3377, only the most extreme anisotropy 
shown (J3 Z ~ 0.3) is able to reduce the gradient of (vf^) 1 ^ 2 along 
the minor axis to the shallow level of the observations. 

To quantify more precisely the findings of the previous para- 
graph, in Fig.[3]we show the plots of A\ 2 = x 2 ~ Xmin> describing 
the agreement between the data and the models, as a function of /3 Z . 
We also show, as a filled square with error bars, the value of the 
global anisotropy we meas ured using more gene ral axisymmetric 
Schwarzschild's models in lCappellari et alj d2007l) . The A\ 2 plots 
confirm the visual impression from the maps of Fig.|2]and the fact 
that the Jeans models give a tight constrain to the anisotropy. The 
small size of the Jeans formal errors is of course illusory, as it is 
due to the restrictive model assumptions. However the comparison 
with the more general results show that the simple anisotropic Jeans 
models are able to recover the anisotropy within the errors. There 
seems to be a systematic difference of A/3 Z ~ 0.05 between the 
anisotropic Jeans and the Schwarzschild determinations, with the 
Jeans results being lower. This difference is at the l evel of the errors 
in the Schwarzschild models (see section 4.3 in ICappellari et al] 
2007), and not particularly surprising given the radically different 
modeling approaches adopted. If real, it may be due to the fact 
that the Jeans models force a constant anisotropy everywhere in 
the galaxy model, while the Schwarzschild ones allow for a more 
realistic and more spherical velocity ellipsoid at intermediate lati- 
tudes (6 ~ 45; see Fig. [TJ. The Jeans models may need a lower 
anisotropy to compensate the mismatch with the data at these inter- 
mediate latitudes. 



4.1.2 M I L comparison 

When studying the stellar population or the dark-matter content of 
galaxies, one is interested more in the M/L determination than 
in the anisotropy. It is important to test whether an error in the 
anisotropy can lead to a significant error in M/L. For this in 
Fig. [4] we show the variation in the dynamical M/L (/-band) for 
each of the models shown in Fig. [2] Also shown for reference 
is the (M/L)schw de termined from the same SAURON data in 
Capp ellari et aD {2006), using three integral Schwarzschild mod- 
els. The comparison shows that the MjL varies very little for the 
ranges of anisotropy observed in fast rotators early-type galaxies. 
In all cases the Jeans M/L at the best-fitting /3 Z represents an im- 
provement over the M / L determined from semi -isotropic models 
(already given in Table 1 o f lCappellari et al.l2006l) . With the excep- 
tion of NGC 5845, the improvement is generally at the level of the 
measurement errors. The latter galaxy is an interesting case on its 



1 The difference between the isotropic models of NGC 821 and 
NGC 3377, and the observations, has already been shown in figure Al of 
ICappellari et~ai] 120071) . 




-0.1 0.0 0.1 0.2 0.3 0.4 

Figure 4. M/L versus anisotropy. The different symbols, connected by 
lines, show the variation with anisotropy f3 z in the dynamical M/L (I- 
band) for each of the Jeans models shown in Fig. [2] Also shown for 
reference, with horizontal lines using the sa me style and color, is the 
(Af/L)g c i lw , determined from the same data in ICappellari etsdl (2006), us- 
ing three-integral Schwarzschild models. The symbols represent NGC 821 
(green circles), NGC 3377 (blue upward triangles), NGC 4621 (black down- 
ward triangles) and NGC 5845 (red squares). The filled symbols correspond 
to the best-fitting f3 z from Fig. fj] The large cross shows the characteristic 
eiTor in /3 Z and (Af/L) Schw . 

own: it has an small size with R c ~ 4'.' 6 (Cappel lari et alj2006h and 
constitute a unique example of an elliptical galaxy with integral- 
field stellar kinematics out to ~ 3R C . Our standard self-consistent 
model still well reproduces the SAURON kinematics without the 
need to invoke a dark matter halo. 

The small dependence of M/L on the anisotropy explains the 
fact that, in a comparison between the semi-isotropic (M / X)j ea ns 
and the three-integral (Af/L)s c hw for a sample of 25 early-type 
galaxies (including the ones in Fig.[4}, we found an excellent agree- 
ment between the two determinations. We did not detect any signif- 
icant bias in the Jeans values and the differe nces could be explained 
by random errors at the 6 per cent rms level ( Cappellar i et al.l2006l . 
their fig. 7). 



4.2 First and second velocity moments model examples 

In the previous section we showed that the anisotropic Jeans for- 
malism introduced in this paper provides a good qualitative de- 
scription of the integral-field second velocity moments observed 
with SAURON. We also showed that for edge-on galaxies the re- 
covered anisotropy agrees well with the one derived from more 
general models. Here we present additional examples of mo del fits 
to the second moments and we show that, using the simple ISatohl 
( 1980) assumption for the splitting of the random and ordered ro- 
tation (Section [3,l,5l l. a remarkably good prediction of the galaxies 
mean velocity can be obtained as well, in many cases. 

In Fig.|5]the data-model comparison is presented for both the 
first and second velocity moments. The model fits were determined 
with the following procedure: 

(i) For each pair of the nonlinear model parameters (i,/3 z ) we 
evaluated the predicted second moments (vf os ) using equa- 
tion <(28j, with b k = 1/(1 - (5 Z ) and Mj = L k .\ 

(ii) The model predictions were scaled by the best-fitting mass- 
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Figure 5. Data-model comparison for three E and three SO fast-rotator galaxies. From left to right, the columns show the fast-rotators NGC 821 (E; i = 
82°, /3 Z = 0.20, K = 0.75), NGC 2549 (SO; i = 90°, /3 Z = 0.17, k = 0.99), NGC 4546 (SO; i = 71°, /3 Z = 0.05, re = 1.00), NGC 4660 (E; 
i = 67°,f3 z = 0.23, re = 0.93), NGC 5308 (SO; i = 87°, f3 z = 0.28, re = 1.02) and NGC 5845 (E; i = 75°, f3 z = 0.18, re = 0.95). The top row 
visualizes the bi-symmetrized SAURON observations for Vims = VV 2 + a 2 . The second row show the best-fitting model predictions for ('"f os ) 1 ^' 2 as given 
by equation 428i . The third row presents the observed SAURON mean velocity V. The bottom row shows the best fitting model first velocity moment v\ os as 
given by equation 438i . The color levels are the same for data and model. These galaxies are all constrained by the photometry to be quite close to edge-on, 
so the models can vary essentially only the single parameter f3 z to fit the shape of the observed first and second moments. The kinematics varies widely for 
different galaxies, yet this single parameter is sufficient to correctly predict the main features of a pair of two-dimensional functions (V r ms and V), once the 
observed surface brightness is given. 



to-light ratio T via equation J5U and the x 2 , describing he agree- 
ment of data and models was determined; 

(iii) Step (i)-(ii) were repeated for a grid of (i,/3 z ) parameters 
and the best fitting value was found. In the grid, the inclination was 
sampled at equally spaced intervals in the intrinsic axial ratio of the 
flattest MGE component, while f3 z was sampled linearly; 

(iv) The Gaussians describing the mass of the MGE (equa- 
tion 1151 ) were scaled by the best fitting T and a prediction for the 
first velocity moment vi^ was computed with equation l l38l l, with 
Kfe = 1; 

(v) The model velocity field «JJJ was scaled by the factor 

— ^n = l Fn\ X nYn\ 

where P is the number of Voronoi bins in the observed SAURON 
data and F n is the corresponding flux. V n and [wu^] n are the ob- 
served and model velocities respectively, and x' n are the bin coor- 
dinates (the x' axis being along the projected major axis). 

The last step ensures that the model has the same projected an- 
gular momentum as the observed galaxy, within the observed re- 
gion. Doing a normal least-squares fit would heavily underesti- 
mate the amount of rotation of the model with respect to the data, 
in cases where the galaxy contains counter-rotating stellar com- 
ponents, which are clearly excluded by the simple assumption of 
a constant k factor for the whole MGE model. The angular mo- 
mentum of the model could be computed by numerically integrat- 
ing L z = vRv^ over the galaxy volume using equation J2U . 
More useful may be to use the fitted galaxy inclination to de- 
projec ted the observed gala xy stellar angular momentum per unit 
mass l lEmsellem et alj [2007). To model kinematically distinct stel- 
lar components one could allow for different components in 



equation l |38l >, and perform a standard least-squares fit as in equa- 
tion d5Tl ). 

For the examples of Fig. [5] we selected three ellipticals and 
three lenticulars fast-rotator showing a range of kinematical prop- 
erties. For all the galaxies the MGE model^j w ere fitted to the ob - 
served photometry using the public software o f lCappellarill2003) . 
Although some differences between the data and models are vis- 
ible, the x 2 P er degree-of-freedom in the fit presented, computed 
from original non- symmetrized data, is close to one for all the fits. 
This indicates that the models are generally consistent with the data, 
within the measurements errors, which is remarkable for what is es- 
sentially a one parameter model! A more general statistical investi- 
gation of the kinematic parameters for a larger sample of galaxies 
goes beyond the scope of the present paper and will be presented 
elsewhere. Here we note that for the galaxies presented in Fig. [5] 
and [6] the parameter k w 1 within a few percent accuracy (ex- 
cept for NGC 821). This appears a general characteristics also of 
the other fast-rotators we modeled. It implies w v\ and con- 
firms the results of more general models. It indicates that the ve- 
locity ellipsoid of fast- rotator early-type galaxies tends to be oblate 
(Cap pellarietal.ll2007l , their fig. 2). 

Interestingly Fig.[5]shows that even the complex double-hump 
structure observed in the velocity fields of NGC 2549, NGC 4660 or 
NGC 5845, which seems to imply a complex dynamical structure in 
the galaxy, can be well reproduced by these one-parameter models, 
once an accurate MGE parametrization for the surface brightness is 
given. The structures are explained by the presence of thin stellar 
disks, already visible in the photometry as strongly disky isophotes. 

15 The MGE parameters for NGC 2549, NGC 4546 and NGC 5308 are 
taken from Scott et al. (in prep aration ) , whil e the MGE models for the other 
galaxies are given in Cappellari et al. (2006). 
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There is no evidence for a significantly different anisotropy between 
the bulge and the disk of these galaxies, as a constant anisotropy 
well reproduces the observations (within \R C ). These examples 
show that a great deal of information on the kinematical structure 
of the fast-rotators is contained in the photometry alone ! 

4.3 Recovery of the galaxy inclination 

We have shown in the previous sections that the anisotropic Jeans 
axisymmetric models presented in this paper, varying essentially 
only one free parameter Z , provide a good qualitative descrip- 
tion of the observed shape of the first and second velocity mo- 
ments of fast-rotator early-type galaxies, when they are already con- 
strained by the photometry to be close to edge-on. We also con- 
firmed that the fast-rotators are characterized by a velocity ellipsoid 
which is flattened in th e z direction so that the anisotropy f3 z > 
( Cappellari et al. I l2007h . In this section we study the variations of 
the models when the inclination is allowed to vary. 

4.3.1 The inclination degeneracy 

In lKrainovic et ai] d2005l) we showed, using real data and analytic 
tests, that there appears to be a degeneracy in the recovery of the 
galaxy inclination (or its corresponding shape) using general three- 
integral axisymmetric models. This is true even in the case of state- 
of-the-art integral-field data, and assuming the galaxy potential is 
accurately known at every inclination and can be uniquely recov- 
ered from the surface brightn ess. This result was co nfirmed for 
a larger sample of galaxies in Cappe llari et al. and in the 

oblate limit, but using a triaxial modeling code, by ?. In a realistic 
situation both the galaxy potential and the stellar luminosity density 
are only approximately known, due to the important intrinsic degen- 
eracy in the deprojection of the surface brightness, due to the possi- 
ble presence of bars, and due to the likely contribution of a (small) 
fraction of dark-matter. Moreover the kinematical data are always 
affected by low-level systematics, which are difficult to control and 
quantify. The dynamical models are also affected by numerical ap- 
proximations and discretisation effects. All this makes it unlikely 
that any inclination derived from dynamical models of the stellar 
kinematics can be trusted, unless further assumptions are made. In 
this section we explore whether the anisotropic Jeans models, com- 
bined with realistic assumptions on the galaxy anisotropy, can be 
used to recover the galaxy inclination. 

4.3.2 Inclination recovery with the Jeans models 

For our experiments on the recovery of the inclination we selected 
the same four galaxies with SAURON integral-fiel d kinematics that 
we pre sented for the same reason in appendix A of Cappellari et al. 
( 2006). These galaxies have a photometry which, under the axisym- 
metric assumption, can be deprojected for a wide range of inclina- 
tions. Most importantly they possess disks of gas or dust, which 
allow the inclination to be determined independently from the stel- 
lar dynamics. In that paper we found that the inclination derived 
using semi-isotropic models appeared in agreement with the one 
from the dust disks. Here we try to understand the reason for the 
good agreement. 

Our prototypical test case is the fast-rotator E4 galaxy 
NGC 2974. The inclination of this galaxy was studied extensively 
by a number of aut hors using the kinematics of the gas (e.g. 
IWeiimans et alj|200 &). From accurate models of the SAURON gas 



kinematics in the reg ions where we have our stellar kinematics 
iKrainovic et al.l {2005) find a best fitting inclination i = 60 ± 3. 
A sequence of model predictions for the velocity second moments 
( u ios) 1,/2 as a function of anisotropy (3 Z and inclination i is shown 
in Fig. [6] A strong similarity is evident between the effect on the 
models of the anisotropy and of the inclination. Another general 
feature one can see from the figure is that the second moments are 
more weakly sensitive to an anisotropy variation when the galaxy 
is intrinsically flat (i — 55°, q ~ 0.32) than when it is rounder 
(i = 90°, q w 0.63). A thin disk can be more strongly anisotropic 
without producing the vertical elongation of the second moments 
(see i = 90°, j3 z — 0.3), which is never observed in real galaxies. 

The observed behavior can be qualitatively understood as fol- 
lows: A more face-on view makes the model intrinsically flatter and 
disk-like, causing the second moments to approach the | cos<^| be- 
havior along the galaxy isophotes as one would expect for a thin 
rotating disk, with a peak on the projected major axis and a strong 
minimum on the minor axis. A similar increase of the second mo- 
ments on the projected major axis, with respect to the minor axis, 
can be produced by reducing the fraction of stars on radial orbits 
(lowering @ z ), thus correspondingly increasing the stars on tangen- 
tial orbits. Due to projection these in fact produce zero motion along 
the minor axis and show a peak on the major axis. For example in 
Fig.[6]the model predictions with (i, f3 z ) = (63.8°, —0.3) are very 
similar to the ones for (i, f3 z ) — (57.5° , 0.0). The contours of the 
X 2 , describing the agreement between the Jeans models and the ob- 
served SAURON V m s are shown in Fig. [7] Within a 3<r confidence 
level any inclination between i — 56° and 61°, and correspond- 
ingly any anisotropy between f3 z = —0.2 and 0.1 is equally consis- 
tent with the data. 

In ICappellari et 

aD fe007l, their fig. 2) we found that all 
the fast-rotators have anisotropy j3 z > 0.05. For the nearly edge- 
on galaxies we confirmed this determination using Jeans mod- 
els in Section 14.1.1 1 and Section 14/21 Forcing this observationally- 
motivated constraint, the inclination of NGC 2974 becomes tightly 
constrained to the range i = 56° to 58°. The qualitative behavior 
seen in NGC 2974 is representative of what we observed for the 48 
galaxies we modeled using this method from the SAURON sample 
(Scott et al. in preparation). In Fig.|7]we show the \ 2 contours for 
the three additional galaxies. In the case of NGC 4150 the incli- 
nation is well constrained by the data independently of anisotropy, 
while in the case of NGC 4459 and NGC 524 the best fitting lo- 
cation has f3 z < and the inclination is heavily dependent on the 
assumed anisotropy. For these two galaxies the un-physical location 
of the x 2 minimum is likely due to the degeneracy in the deprojec- 
tion of the surface brightness (Section [3.1.21 ). These two galaxies 
have a quite low inclination and the deprojected mass distribution is 
probably significantly in error even at the true inclination, which is 
then formally excluded by the data. Like all the fast-rotators, these 
galaxies likely possess stellar disks, which cannot be detec ted at 
these low inclinations teix & Whiteil 199d : lMagorriar] [r999). This 
would explain the general fact that we tend to better reproduce the 
kinematics of the fast-rotators which are close to edge-on than the 
kinematics of the more face-on ones. Nonetheless, in all the four 
cases, if one enforces the constrain (3 Z > 0.05, the recovered incli- 
nation lies within the error bars of the independent gas/dust deter- 
minations. This indicates that one can use these anisotropic Jeans 
models to recover the inclinations of the fast-rotators by enforc- 
ing the requirement f3 z > 0.05, and this explains the good agree- 
ment with the determi nations from isotropic models we found in 
ICappellari et alj (2006). 
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Figure 6. Inclination and anisotropy variation for NGC 2974. Each of the five columns show the model predictions for the second velocity moment {vf^) 1 ^ 2 
computed at different values of f} M = —0.3, —0.15, . . . , 0.3. Different rows correspond to different inclinations i, equally spaced in the axial ratio of the 
flattest Gaussian in the MGE model of NGC 2974. The model mass-to-light ratio T was optimized for each fit using equation 15 1 1 - The values for the input 
parameters (/3 Z , i) and for the best fitting T are printed on top of each panel. The bottom row shows (i) the bi-symmetrized S AURON data V rms = VV 2 + cr 2 ; 
(ii) the best fitting model for (vf^) 1 / 2 ', (hi) the bi-symmetrized SAURON data for V and (iv) the best-fitting model for v\^. The colormap has the same scaling 
for all the panels. 



4.4 Relation with the tensor virial theorem 

The Jeans equations are useful when one needs a predictions for the 
spatial variation of the kinematics in galaxies. However the stars 
in galaxies also satisfy more general global relations between kine- 
matic observables, the most imp ortant of which are the virial equa- 
tions jBinnev & Tremainell 19871 . sec. 4.3). These are a tensor gen- 
eralization of the well known virial theorem 2K — —W, which 
relates the total kinetic energy K and the potential energy W of an 
isolated stellar system in equilibrium. 

The small dependence of the measured Mf L on the anisotropy 
(Section ^. i.2\ can be qualitatively understood as a consequence of 
the tensor virial equations. For an axisymmetric galaxy in equilib- 
rium the total kinetic energy is proportional to the galaxy mass. The 
ratio of the kinetic energy K zz along the direction of the symme- 
try z-axis and the kinetic energy K xx along a direction orthogo- 
nal to it, is only a function of the galaxy shape (equation [4-91] of 
iBinnev & Tremainelll987l) . The effect of increasing the anisotropy 
f3 z is simply that of redistributing the same kinetic energy K xx 
from the tangential to the radial direction. This reduces the veloc- 



ity second moments along the galaxy projected major axis, while 
correspondingly increasing them along the minor axis. The pro- 
jected luminosity-weighted second moments, from which the M/L 
is measured, remain nearly unchanged when averaged over the full 
galaxy image. 

The virial equations also provide a qualitative understanding of 
why the observations can constrain the inclination of an axisymmet- 
ric galaxy (Section[43]2). In fact they state that the ratio K xx /K zz 
increases with galaxy flattening. In other words, for a given pro- 
jected surface brightness, a flat system needs more kinetic energy 
in a direction parallel to its the equatorial plane, to be in hydrostatic 
equilibrium. 

The fast rotator NGC 2974 of Fig. [6] shows high velocity sec- 
ond moments along the projected major axis. Although the galaxy 
does not appear very flat in projection, its kinematics is already an 
indication that the system is likely to be intrinsically flat. A flat 
model has more kinetic energy in the equatorial plane and will be 
able to easily reproduce the high V rms along the major axis, with- 
out invoking extreme tangential anisotropy. A quantitative example, 
using both the first and second velocity moments, of how the ten- 
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Figure 7. The inclination-anisotropy degeneracy. Contours of A\ 2 = 
X 2 — Xmin' describing the agreement between the data and the anisotropic 
Jeans models, are plotted as a function of the anisotropy parameter j3 z and 
the galaxy inclination i in degrees. The lowest levels correspond to the for- 
mal 1, 2 and 3tr (thick red line with A\ 2 = 11.8) confidence levels for two 
degrees of freedom. Additional contours are characterized by a factor of 2 
increase in A\ 2 - The two horizontal blue solid lines enclose the region of 
acceptable inclinations as determined from the gas kinematics or dust ge- 
ometry (appendix A o flCa ppellari et al.|[2 0Q6). The vertical solid green line 
indicates isotropy (/3 Z = 0). The different panels correspond to NGC 4150, 
NGC 2974, NGC 4459 and NGC 524 respectively. There is a clear degen- 
eracy between inclination and (3 Z anisotropy. However if one enforces the 
observationally-motivated constraint /3 2 > 0.05, the recovered inclination 
always lies within the error bars of the gas/dust determinations. 



sor virial equations can be used to constrain the galaxy inclination 
is presented in fig. 1 1 o f Capp ellari et"ai1 j2007l) . There we showed 
that the (V/ a, e) diagram of Binnev ( 2005) can be used to provide a 
qualitative estimate of the inclination of fast-rotators galaxies, given 
some assumptions on the anisotropy. 

Although the virial equations can be used to estimate the M/L 
and the inclination of galaxies, they have some obvious limitations: 
(i) They do not allow one to make full use of the spatial informa- 
tion contained in the integral-field observations. The whole two- 
dimensional kinematics is reduced to one number; (ii) They do not 
allow the limited spatial extent of the kinematics to be properly 
taken into account; (iii) They do not easily deal with non-similar 
isophotes and multiple photometric components in galaxies. 



5 SUMMARY 

We present a generalization of the widely used semi-isotropic (two- 
integral) axisymmetric Jeans modeling method to describe the stel- 
lar dynamics of galaxies. Our method uses the powerful Multi- 
Gaussian Expansion (MGE) technique to accurately parameterize 
the galaxies photometry. It represents an anisotro pic extension of 
what w as presented in the semi-isotropic case by Emsellem et al. 
j 1994b . and it maintains its simplicity and computational efficiency. 



We assume (i) a constant mass-to-light ratio and (ii) a velocity el- 
lipsoid which is aligned with the cylindrical (R, z) coordinates and 
has a flattening quantified by the classical z-anisotropy parameter 
ftz = 1 — Vz/vjj, where z is the galaxy symmetry axis. 

We test the technique usin g SAURON integral-fi eld observa- 
tions of the stellar kinematics temsellem et al.l 2004) for a small 
set of fast-rotator galaxies with a variety of kinematical proper- 
ties. For galaxies that are constrained by the photometry to be 
close to edge-on we find that, although in the semi-isotropic limit 
((3 Z = 0) the models do not provide a good description of the data, 
the variation of the single global anisotropy /3 Z is generally suffi- 
cient to accurately predict the shape of both the first (V) and sec- 
ond (Vims = y/V 2 + a 2 ) velocity moments, once an detailed MGE 
parametrization of the photometry is given. An accurate descrip- 
tion of the photometry, including ellipticity variations and disky 
isophotes, appears crucial to reproduce in detail the features of the 
kinematics. 



In all cases we find that /3 Z > 0.05, while generally a 



with good accuracy. This confirms previous findings on the dynam- 
ical structure of these galaxies, showing that their velocity ellip- 
soid tends to be oblate dCappellari et alj 2007t) . The anisotropy de- 
rived with our anisotropic Jeans dynamical modeling method agrees 
within the errors with the one previously measured using a more 
general Schwarzschild approach. 

For fast-rotator galaxies that are not constrained by the pho- 
tometry to be close to edge-on, we find that in general the incli- 
nation i (or the corresponding galaxy shape) and the anisotropy 
j3 z are highly correlated and cannot be independently determined. 
However, if we introduce the observationally-motivated constraint 
f3 z > 0.05, the inclination becomes constrained to a narrow range 
of values and it agrees with independent determinations when those 
are available. 

We are applying this method to determine the inclination, the 
mass-to-light ratio and the amount of rotation of a large sample of 
galaxies for which integral field kinematics are available. For galax- 
ies close to edge-on, the global anisotropy or the dynamical struc- 
ture of different galaxy subcomponents (e.g. bulge and disk) can 
also be investigated. We are using this method to test independent 
determinations of the masses of supermassive black holes in the nu- 
clei of fast-rotator galaxies. This technique is ideal to study the dark 
matter content and the anisotropy of disks of spiral galaxies. 
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erf 



(Ly/2) - y 



+ erf 



(Ly/2) + y 



(A6) 



In practice the integral of equation ( IA5b can be limited to the region 
where K(x, y) is significantly nonzero: 

- (L x /2) - 3ma,x{a z }<x<(L x /2) + 3max{<Ti} (A7) 
-(L y /2) - 3max{ai} <y < (Ly/2) +3max{a,}. (A8) 



APPENDIX A: SEEING AND APERTURE CONVOLUTION 

The first and second velocity moments have to be convolved with 
the instrumental point-spread-function (PSF) and integrated over 
the aperture used in the observations, before making comparisons 
with the observed quantities. The PSF effects are particularly im- 
portant when studying the nuclear regions of galaxies (e.g. to esti- 
mate the mass of supermassive black holes). The PSF-convolved 
velocity moment s are given by the ge neral formulas (e.g. equa- 
tions [51]-[53] of lEmsellem et alii 994) 

E obs = E®PSF (Al) 
r i (ggg) g PSF 

blos]obs = ^ (A2) 

-g- _ (S<)®PSF 

RoJobs - = , (A3) 

^obs 

with Cg> representing convolution and where the normalized PSF is 
generally described as the sum of Q Gaussians 

PSF(*) = E^exp(^), (A4) 

with YliLi Gi = 1. To efficiently evaluate numerically the above 
convolutions up to large radii, and even for very flat models, one can 
evaluate the models predictions Eu 2 os and ETJITjI on a grid linear in 
the logarithm of the elliptical radius and in the eccentric anomaly. 
This is done by defining a logarithmically-spaced radial grid Rj and 
then computing the moments at the coordinate positions (x' , y') — 
(Rj cos 9k, q' Rj sin 9k), for linearly spaced 9k values, with q' a 
characteristic (e.g. the median) observed axial ratio of the MGE 
model. The model is then re interpolated onto a fine Cartesian grid 
before the convolution with the PSF using fast Fourier methods. 
Finally the model is rebinned into the observed apertures. 

Alternatively, especially when the apertures sparsely sample 
the plane of the sky, instead of performing the convolutions of equa- 
tions lAlHA3l on a regular grid, one can evaluate the same convolu- 
tion integrals only for the observed apertures, w hile also includin g 
the integration onto the apertures (appendix D of Oia n et aO 19951) . 
Given a rectangular aperture of sides L x and L y , aligned with the 
coordinate axes at position (x' , y') on the sky, and assuming a PSF 
given by equation l lA4t . the PSF-convolved observable S(x',y') 
inside the aperture is 

r CO p CO 

S b s (x',y')= / S(x',y') K(x - x',y - y')dxdy, (A5) 



where S has to be replaced by E, Eui os , or Euj 2 os respectively, and 
correspondingly Sobs becomes E obs , [E75io7]obs, or [Ew^Jobs re- 
spectively. The kernel is given by 

'-' " ' r " - 1 r (L x /2)+x 



K(x,y) = E^jerf 



(W2) - 
V2a t 



+ erf 



V2a t 
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